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Abstract. The proton radioactivity lifetimes of spherical proton emitters from the ground and the 
isomeric states are calculated using the microscopic nucleon-nucleus interaction potentials. The 
daughter nuclei density distributions are folded with a realistic density dependent M3Y effective 
interaction supplemented by a zero-range pseudo-potential. The density dependence parameters of 
the interaction are extracted from the nuclear matter calculations. The saturation energy per nucleon 
used for nuclear matter calculations is determined from the co-efficient of the volume term of Bethe- 
Weizsacker mass formula which is evaluated by fitting the recent experimental and estimated atomic 
mass excesses from Audi- Wapstra-Thibault atomic mass table by minimizing the mean square devi- 
ation. The quantum mechanical tunneling probability is calculated within the WKB approximation. 
Spherical charge distributions are used for calculating the Coulomb interaction potentials. These 
calculations provide good estimates for the observed proton radioactivity lifetimes. 
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1. Introduction 

The discovery of proton radioactivity has provided a tool to obtain spectroscopic informa- 
tion [1], [2] near the proton drip line. The decaying proton is the unpaired one not filling 
its orbit. These decay rates are sensitive to the proton decay Q values and the orbital an- 
gular momenta which in turn help to determine the orbital angular momenta of the emitted 
protons. It has the lowest Coulomb potential among all possible charged clusters in nuclei 
and mass being the smallest it suffers the highest centrifugal barrier, enabling this process 
suitable to be dealt within WKB barrier penetration model. In the present work we provide 
estimates for the proton radioactivity lifetimes of some spherical proton emitters from the 
ground and the isomeric states using the nucleon-nucleus interaction potentials obtained 
microscopically by single folding the daughter nuclei density distributions with a realistic 
density dependent M3Y effective interaction (DDM3 Y) whose density dependence param- 
eters have been extracted from the nuclear matter calculations. 
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2. Energy coefficients of Bethe-Weizsacker mass formula from atomic mass excesses 

The Bethe-Weizsacker formula [3] is an empirically refined form of the liquid drop model 
for the binding energy of nuclei. Expressed in terms of the mass number A and the atomic 
number Z for a nucleus, the binding energy B(A,Z) is given by, 

^ .on Z(Z-l) (A-2Z) 2 
B(A, Z) = a v A - a s A 2 / 3 - a c V ^ 1/3 ' - a sy J- -j-L- + S, (1) 

where S = a p A~ 1/2 , for even N-even Z, = -a p A~ 1/2 for odd N-odd Z, = for odd A 
and N is the neutron number of the nucleus. Using the definition of nuclear binding energy 
B(A,Z) which is defined as the energy required to separate all the nucleons consituting a 
nucleus, the mass equation can be expressed as 

M nucleus (A, Z) = Zm p + (A - Z)m n - B(A, Z), (2) 

where m p and m n are the masses of proton and neutron respectively and M nuc i eus (A, Z) 
is the actual mass of the nucleus. The mass of the nucleus can be obtained from the atomic 
mass excess by use of the relationship 

M nucleus (A, Z) = Au + AM A , Z - Zm e + a el Z 2m + b e iZ 5 35 , (3) 

where m e is the mass of an electron, the atomic mass unit u is 1/12 the mass of 12 C atom, 
AMa z is the atomic mass excess of an atom of mass number A and atomic number Z 
and the electronic binding energy constants [4] a e j = 1.44381 x 10~ 5 MeV and b e i = 
1.55468 x 1CP 12 MeV. Hence from Eq.(2) and Eq.(3) the atomic mass excess is given by 

AM AiZ = ZAm H + (A - Z)Am n - a el Z 2 39 - b el Z 5 35 - B(A, Z) (4) 

where Am# = m p + m e — u = 7.28897050 + a e i + b e i MeV and Am„ = m„ — u = 
8.07131710 MeV. Root mean square deviation a for atomic mass excesses is defined as 

a 2 = (1/JV) J2l AM Th - AM Ex ] 2 (5) 

where the summation extends to N data points for which experimental atomic mass ex- 
cesses are known. The quantity AMxh is theoretically calculated atomic mass excess 
obtained using Eq.(4) while AM& is the corresponding experimental atomic mass ex- 
cess obtained from the Audi-Wapstra-Thibault atomic mass table [5]. For the calculations 
of the mean square deviation <r 2 , masses [5] of 3179 nuclei including the 951 extrapo- 
lated values which are predicted according to systematics, have been used for least square 
fitting. The values of the energy coefficients are a v — 15.2602 ± 0.019671^6^, a s — 
16.2673±0.0621MeU, a c = 0.688974±0.001357MeU, a sym = 22.2090±0.0481MeV, 
a p = 10.0757 ±0.8536MeV. 



3. The volume energy coefficient and the nuclear matter incompressibility 

A density dependent M3Y effective nucleon-nucleon (NN) interaction [6] based on the G- 
matrix elements of the Reid-Elliott NN potential has been used to determine the nuclear 
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matter equation of state. The equilibrium density of the nuclear matter has been determined 
by minimizing the energy per nucleon. The density dependence parameters have been 
chosen to reproduce the saturation energy per nucleon and the saturation density of spin 
and isospin symmetric cold infinite nuclear matter, which is called the standard nuclear 
matter. The general expression for the DDM3Y interaction potential v(s) is written as 

v(s,p,e) = t M3Y (s,e)g(p,e) = Ct M3Y (1 - (3(e)P V3 ) (6) 

where M3Y effective interaction potential supplemented by a zero range pseudopotential 
t M3Y is given by [7] 

—4s —2.5s 

;M3Y = 799 9 f_ _ 2134^- + Joo(E)5(s) (7) 

where the zero-range pseudo-potential Joo(e) representing the single-nucleon exchange 
term is given by 

Joo(e) = -276(1 - ae)(MeV.fm 3 ) (8) 

The energy per nucleon e obtained using the effective nucleon-nucleon interaction v(s) for 
the standard nuclear matter, is given by 

= 3H 2 k% g(p,e) P J v = 3ft 2 fc 2 P J v C(l - Pp 2 ' 3 ) 

1 10m J 2 [ 10m J 1 2 J K ' 

where m is the nucleonic mass, k F = (1.5tt 2 p) 1 / 3 is the Fermi momentum, p is the nucle- 
onic density while po being the saturation density for the standard nuclear matter and J v 
represents the volume integral of t M3Y , the M3Y interaction supplemented by the zero- 
range pseudopotential. Eq.(9) can be differentiated with respect to p to yield equation 

de/dp = [fi 2 k 2 F /hmp] + J V C[1 - (5/3)/3p 2/3 ]/2 (10) 

The equilibrium density of the nuclear matter is determined from the saturation condition 
de/dp = 0. Then Eq.(9) and Eq.(10) with the saturation condition can be solved simulta- 
neously for fixed values of the saturation energy per nucleon e and the saturation density 
po of the standard nuclear matter, to obtain the values of the density dependence parameters 

(5= [(l-p)p- 2/3 ]/[3-(5/3)p], p= [10 meo ]/[h 2 k 2 Fo ], (11) 



C = -[2n 2 k 2 Fo ]/[5mJ vPo (l - (5/3)/3po /3 )] with k F „ = [1.57r 2 p ] 1/3 , (12) 
respectively. The incompressibility K n of the standard nuclear matter is given by 

K = k 2 F d 2 e/8k 2 F = 9p 2 8 2 e/dp 2 | p=Po = [-(^f 1 ) - 5J v C^p 5 /3 ] (13) 

Identifying a v as the saturation energy per nucleon for the spin and the isospin symmet- 
ric cold infinite nuclear matter and hence using the saturation energy per nucleon equal to 
— 15.2602MeV r along with the density dependent M3Y effective interaction and the com- 
monly used values for saturation density equal to 0.1533/m~ 3 [8] and energy dependence 
parameter a = O.OOHMeV^ 1 , the density dependence parameters have been found to be 
C = 2.07, (3 = 1.668/m 2 while the nuclear incompressibility is found to be 293.4MeV r 
which is in close agreement with experimental data [9] and theoretical estimates [8]. 
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4. The half lives of the spherical proton emitters 

The microscopic nuclear potentials Vn [R) have been obtained by single folding the den- 
sity of the daughter nucleus with the finite range realistic DDM3Y effective interacion 
as 

V N (R) = J p(r)v[\r-R\]d 3 r (14) 

where R and f are, respectively, the co-ordinates of the emitted proton and a nucleon be- 
longing to the residual daughter nucleus with respect to its centre. The density distribution 
function p used for the daughter nucleus, has been chosen to be spherically symmetric 

p(r) = p /[l + exp({r - c)/a)] (15) 



where c = r p (l - TT 2 a 2 /3r 2 p ), r p = 1A3A 1 / 3 , a = 0.54 fm (16) 

and the value of p a is fixed by equating the volume integral of the density distribution 
function to the mass number Ad of the residual daughter nucleus. The distance s between 
any nucleon belonging to the residual daughter nucleus and the emitted proton is given by 
s = \ f— R\ while the interaction potential between these two nucleons v(s) appearing in 
Eq.(14) is given by the DDM3Y effective interaction. The total interaction energy E(R) 
between the proton and the residual daughter nucleus is equal to the sum of the nuclear 
interaction energy, Coulomb interaction energy and the centrifugal barrier. Thus 

E(R) = V N (R) + V C (R) + h 2 l(l + l)/(2pR 2 ) (17) 

where p = nipirid/mA is the reduced mass, m p , irid and itia are the masses of the pro- 
ton, the daughter nucleus and the parent nucleus respectively, all measured in the units of 
MeV/c 2 . Assuming spherical charge distribution for the residual daughter nucleus, the 
proton-nucleus Coulomb interaction potential Vc (R) is given by 

Vc{R) = {^-)\i~{^) 2 ]forR<R Cl =^- otherwise (18) 
ZJx c Jx c Jx 

where Zd is the atomic numbers of the daughter nucleus. The touching radial separation R c 
between the proton and the daughter nucleus is given by R c = c p +Cd where c p and Cd have 
been obtained using Eq.(16). The energetics allow spontaneous emission of protons only if 
the released energy Q is a positive quantity. The half life of a parent nucleus decaying via 
proton emission is calculated using the WKB barrier penetration probability. The assault 
frequency v is obtained from the zero point vibration energy E v = (l/2)hv. The decay 
half life T of the parent nucleus (A, Z) into a proton and a daughter (Ad, Zd) is given by 

r=[(ftIn2)/(2£7„)][l + exp(A')] (19) 

where the action integral K within the WKB approximation is given by 

K = (2/h) / [2p(E(R) -E v - Q)] 1/2 dR (20) 
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where R a and Rt, are the two turning points of the WKB action integral determined from 
the equations 

E(R a ) = Q + E v =E(R b ). (21) 

From a fit to the experimental data on cluster emitters a law given by eqn.(5) of ref. [10], 
which relates E v with Q, was found. For the present calculations the same law [10] ex- 
tended to protons for the zero point vibration energies along with the experimental Q values 
[1], [2] have been used. The same set of data of ref. [2] has been used for the present cal- 
culations and presented in table- 1 below. In fig. 1, the plots of theoretical and experimental 
proton radioactivity half lives have been presented. 

Table 1. Comparison between experimentally measured and theoretically calculated 
half-lives of spherical proton emitters. The asterisk symbol (*) denotes the isomeric 
state. The experimental Q values, half lives and I values are taken from reference [1], 
The results of the present calculations have been compared with the experimental values 
and with the results of UFM [1] and DWBA [2] estimates. 
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Ang. 


Released 


Experimental 


Present 


UrM 


nucleus 


mom. 


Energy 


logio |hail Irte J 


calculations 


DWrSA 


Symb. 


i(h) 


Q(MeV) 


logwT(s) 


log w T(s) 


log w T(s) 


109 j 


2 


0.829 


-4.000 


-4.25 


-5.000 ft 


112 Cs 


2 


0.823 


-3.301 


-3.13 


-4.167 ft 


113 Cs 


2 


0.977 


-4.770 


-5.53 


-6.26% b 


146 Tm 


5 


1.140 


-0.629 


0.52 


-0.456 6 


i46 Tm * 


5 


1.210 


-1.143 


-0.30 


-l.216 b 


147 Tm 


5 


1.071 


0.591 


0.99 


1.095" 


147 Tm* 


2 


1.139 


-3.444 


-3.38 


-3.199" 


150 Lu 


5 


1.283 


-1.180 


-0.58 


-0.859" 


161 Lu 


5 


1.255 


-0.896 


-0.65 


-0.573" 


156 Ta 


2 


1.028 


-0.620 


-0.38 


-0.461" 


156 Ta . 


5 


1.130 


0.949 


1.67 


1.446" 


157 Ta 





0.947 


-0.523 


-0.43 


-0.126" 


160 Re 


2 


1.284 


-3.046 


-2.99 


-3.109" 


161 Re 





1.214 


-3.432 


-3.46 


-3.231" 


161 Re* 


5 


1.338 


-0.488 


-0.58 


-0.458" 


165 Jr * 


5 


1.733 


-3.469 


-3.51 


-3.428" 


166 /r 


2 


1.168 


-0.824 


-1.10 


-1.160" 


166 j r » 


5 


1.340 


-0.076 


0.23 


0.021" 


167 Jr 





1.086 


-0.959 


-1.27 


-0.943" 


167 Jr * 


5 


1.261 


0.875 


0.71 


0.890" 


m Au* 


5 


1.718 


-2.654 


-3.03 


-2.917" 


1S5 Bi 





1.624 


-4.229 


-5.43 


-5.184" 


185 Bi* 


5 


1.611? 


-4.357 


-1.43 


-1.678 6 
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Figure 1. Comparision between theoretical and experimental proton radioactivity half 
lives. Dots represent the observed half lives while the circles and the triangles represent 
those corresponding to the present and the UFM[1]/DWBA[2] calculations respectively. 

5. Summary and conclusion 

The half lives for proton-radioactivity have been analyzed with microscopic nuclear poten- 
tials obtained by the single folding the DDM3Y effective interaction. The energy depen- 
dence parameters of the interaction have been obtained from the nuclear matter calcula- 
tions. The saturation energy per nucleon used for nuclear matter calculations is determined 
by fitting the theoretical mass excesses based on semi-empirical liquid drop model to the 
experimental mass excesses. This procedure of obtaining nuclear interaction potentials 
has a profound theoretical basis. The results of the present calculations for the proton- 
radioactivity lifetimes are in good agreement over a wide range of experimental data. 
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